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ABSTRACT 

Context. With the increasing knowledge of the terrestrial planets due to recent space probes it is possible to model their rotation 
with increasing accuracy. Despite that fact, an accurate determination of Venus precession and nutation is lacking 
Aims. Although Venus rotation has been studied in several aspects, a full and precise analytical model of its precession-nutation 
motion remains to be constructed. We propose to determine this motion with up-to-date physical parameters of the planet 
Methods. We adopt a theoritical framework already used for a precise precession-nutation model of the Earth, based on a 
Hamiltonian formulation, canonical equations and an accurate development of the perturbing function due to the Sun. 
Results. After integrating the disturbing function and applying the canonical equations, we can evaluate the precession constant 
^ and the coefficients of nutation, both in longitude and in obliquity. We get 'l/ = 4474" .35/ Jcj/ ± 66.5, corresponding to a 
precession period of 28965.10 ± 437 years. This result, based on recent estimations of the Venus moment of inertia is signifi- 
cantly different from previous estimations. The largest nutation coefficient in longitude with an argument 2Ls (where Lg is the 
longitude of the Sun) has a 2" 19 amplitude and a 112.35 d period. We show that the coefficients of nutation of Venus due to 
its triaxiality are of the same order of amplitude as these values due to its dynamical flattening, unlike of the Earth, for which 
they are negligible. 

Conclusions. We have constucted a complete theory of the rotation of a rigid body applied to Venus, with up-to-date determi- 
nations of its physical and rotational parameters. This allowed us to set up a new and better constrained value of the Venus 
precession constant and to calculate its nutation coefficients for the first time. 

Key words. Venus rotation 



1. Introduction 

Among the planets of the Solar system, Venus shows pe- 
culiar characteristics : its rotation is retrograde and very 
slow. Since its 243.02 d period was determined by radar 
measurements (Golstein 1964 , Carpenter 1964), several 
authors have attempted to explain these two character- 
istics. Goldreich and Peale (1970) showed that thermally 
driven atmospheric tidal torques and energy dissipation at 
the boundary between a rigid mantle and a differentially 
rotating liquid core are possible mechanisms to maintaine 
the retrograd spin. Lago and Cazenave (1979) studied the 
past evolution of the rotation of Venus using the hypothe- 
sis that only solar tidal torques and core-mantle coupling 
have been active since its formation. They found it con- 
ceivable that Venus originally had a rotation similar to the 
other planets and evolved for 4.5 x 10^ years from a rapid 
and direct rotation to the present slow retrograde one. 
Others authors proposed different scenarios by suppos- 
ing a high value of the initial obliquity (Dobrovoskis,1980, 
Mc Cue and Dormand,1993). Yoder (1995) gave a full 
account of the various internal mechanisms acting on 
Venus'obliquity, such as core friction, CMB (core-mantle 
boundary) ellipticity, and resonant excitations. Correia 



and Laskar (2001,2003 and Correia et al., 2003), explored 
a large variety of initial conditions in order to cover possi- 
ble formation and evolutionary scenarios. They confirmed 
that despite the variations in the models, only three of the 
four final spin states of Venus are possible and that the 
present observed retrograd spin state can be attained by 
two different processes : one prograde and the other ret- 
rograde. Although these various studies concern the evo- 
lution of Venus rotation on very long times scales, few 
attempts have been made to accurately model the rota- 
tion for short times scales. An analytical study of the ro- 
tation of a rigid Venus model carried out by HabibuUin 
(1995) with rather uncommun parametrization, showing 
that Venus rotates almost uniformly with negligible libra- 
tion harmonics. In the following, we propose an alternative 
construction of a rigid Venus rotation model, contradicting 
to some extent these last results. First we show the char- 
acteristics of Venus with respect to the corresponding ones 
for the Earth (section [2]) . We explain the parametrization 
of the rotation of a rigid Venus using the Andoyer variables 
(section ^ . We then determine the reference points and 
planes to study the rotational dynamics of Venus (section 
E]). We give the equations of motion and the analytical 
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developments that allow us to obtain the precession and 
the nutation of the planet (sections [5l|6] and [7]). Finally we 
give the precession, the tables of nutation and the polar 
motion of Venus (section |H] and We will discuss our 
results in section (TU] In this study we make some approx- 
imations : we suppose that the relative angular distance 
between the three poles (the poles of angular momentum, 
of figure and of rotation) is very small as is the case for the 
Earth. Moreover we solve a simplified system where the 
second order of the potential does not appear. Last, we 
suppose that the small difference between the mean longi- 
tude and the mean anomaly of the Sun which corresponds 
to the slow motion of the perihelion, can be considered as 
constant. In sections [3] to [TUl the orbital elements of Venus 
and their time dependence are required to carry out our 
numerical calculations. In these sections we used the mean 
orbital elements given by Simon et al.(1994). These ele- 
ments are valid over 3000 years with a relative accuracy 
of 10~^. This prevents our theory being used in long term 
calculations, our domain of validity being 3000 years. 

2. Venus characteristics 

From the point of view of its physical characteristics, 
Venus can be considered as similar to sister of the Earth: 
the two planets have roughly the same size, mass and mean 
density (see Table [TJ). 





Earth 


Venus 


Semi-major axis 


1.000001 U.A 


0.723332 U.A 


Eccentricity 


0.0167 


0.0068 


Inclination 





3.39 


Ascending node 


174.87 


7667 


Period of revolution 


365.2,5 d 


224.70 d 


Density 


5.51 g/cm^ 


5.25 g/cm^ 


Equatorial radius 


6378 km 


6051 km 


Mass 


5.97 * 10" kg 


4.87 * 10='"' kg 


Obliquity 


23.43 


2.63 


Period of rotation 


0.997 d 


-243.02 d 


Triaxiality : 


-5.34* 10"" 


-1.66 t 10"" 


Dyn. flattening : ^°2C~"° 


3.27*10"-' 


1.31 * lO"-"^ 


B-A 
2[2C-A-B) 


0.0016 


0.12 



Table 1. Comparison between the orbital (at J2000.0) 
and physical parameters of Venus necessary for our study 
and the corresponding ones for the Earth. 



Despite this fact, the rotation of the two planets has 
with h distinctive characteristics. First, the rotation of 



Venus is retrograde with a period of 243 days. Furthermore 
the obliquity of Venus is very small (2.63) compared to 
that of the Earth (23.44). Last, as we will see in the fol- 
lowing, the perturbating function for the rotational mo- 
tion of the planets depends on the dynamical flattening 
Hv = ^'^~2c~^ and on the triaxiality Ty = The 
Earth has a fast rotation, therefore it has a relatively 
large dynamical flattening. Moreover its main moments 
of inertia A and B according to the X and Y axes can be 
considered equal when comparing them to C . Thus the co- 
efficient of triaxiality is very small with respect to the coef- 
ficient of dynamical flattening. On the other hand, Venus 
has slow rotation. Therefore its coefficient of dynamical 
fiattening is less important than that of the Earth but 
its coefficient of triaxility is of the same order as the co- 
efficient of dynamical flattening (see Tabled]). For these 
reasons, it is important to account for the effect of triax- 
iality on the rotation of Venus and it is interesting to see 
what difference will emerge with respect to the Earth's 
rotation. In Table [TJ for the sake of clarity we have writ- 
ten the orbital parameters at J2000.0. However through- 
out this paper the orbital parameters are considered as a 
function of time (containing terms of first and third order 
in time). 

3. Description of the motion of rotation and 
torque-free motion of Venus 

In order to describe in the most complete way the rotation 
for a rigid model of Venus around its center of mass we 
need four axes : an inertial axis (arbitrairly chosen) , the 
figure axis, the angular momentum axis, and the axis of 
rotation (see FiglTJ). 




Fig. 1. Relation between the Euler angles and the 
Andoyer variables. 

• As an inertial plane (0, X, Y) we choose the Venus 
osculating orbital plane at the epoch J 2000.0. We note 
Toy, the intersection between this orbital plane and the 
mean equator of Venus at J2000.0, which is the reference 
point on the (0, X, Y) plane. The choice of this plane is 
natural because we want to compute the variation of the 
obliquity of Venus which is the angle between the orbit of 
Venus and the plane normal to L. 
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• The angular momentum axis of Venus is tfie axis 
directed along (f , with G the amplitude of the angular 
momentum. We denote h — TgyQ and I as respectively 
the longitude of the node with respect to Tgy and the 
inclination with respect to (0,X,Y). g and J are the longi- 
tude of the node and the obliquity of this same plane with 
respect to Q and to the equatorial plane . 

• The figure axis is the principal axis of Venus directed 
along Oz and perpendicular to the equatorial plane. We 
choose it so that it coincides with the axis of the largest 
moment of inertia C . We denote Ox the axis that is pointed 
toward the origin meridian on Venus, which can be defined 
in a conventional way. The axis of the figure is determined 
with respect to the inertial plane (0,X,Y) through the 
Euler angles hf,If. The parameter $ « l + g + h gives the 
position of the prime meridian (0,x) with respect to Toy. 

• The axis of rotation is determined through the vari- 
ables hr,Ir, i-C- the longitude of the node and the inclina- 
tion with respect to Tpy and the orbit of Venus at J2000.0. 
For the sake of clarity, these parameters are not repre- 
sented in FigHJ 

To describe the motion of the rotation of the 
rigid Earth, Kinoshita (1977) used Andoyer variables. 
Therefore we have introduced these variables to apply an 
analogous theory to Venus. The Andoyer variables (1923) 
are (see FigH]) : 

— L the component of the angular momentum along the 
Oz axis 

— H the component of the angular momentum along the 
OZ axis 

— G the amplitude of the angular momentum of Venus 

— h the angle between Toy and the node Q between 
the orbital plane and the plane normal to the angular 
momentum. 

— g the angle between the node Q and the node P 

— / the angle between the origin meridian Ox and the 
node P. 

This yields: 

L^G-cosJ et H ^G-cosI (1) 

where /, J are respectively the angle between the angular 
momentum axis and the inertial axis (0,Z), between the 
angular momentum axis and the figure axis. 

By using the spherical trigonometry in the triangle (P, 
Q, R) (see FigHJ), we determine the relations between the 
variables : 

cos// = cos /cos J — sin /sin J cos g (2) 

sm{hf — h) sin($ — /) sing 
sin J sin/ sin// 

The angle J is supposed to be very small, so by neglecting 
the second order we obtain (Kinoshita, 1977): 

hf = h+ -^smg + 0{J^) (4) 



If = I + Jcosg + 0{J^) (5) 

$ = / + .g- Jcot/sin5 + 0(j2). (6) 

hf and // correspond to the same definition as h and /, 
but for the axis of figure instead of the axis of angular 
momentum. The Hamiltonian for the torque-free motion 
of Venus corresponding to the kinetic energy is : 

^ l.sin^Z cos^Z,,^2 t2^ 1-^^ 
^°=2(^ + ^)(^ -^)+2C 

where A, B, G are the principal moments of inertia of 
Venus. Then the Hamiltonian equations are: 

|(L,G,//) = -^|^ 

d , - dFo 

dt^^^'^'^^^di^Gjry 

The Hamitonian is the full energy of the system. We see 
the Hamiltonian ([7]) corresponds to the kinetic energy. 

4. Hamiltonian referred to a moving plane 

This section is based on the pioneering work by Kinoshita 
(1977) for the rotation of the Earth. As a new reference 
plane, we choose the plane of the orbit of Venus at the date 
t instead of the fixed orbital plane at J2000.0. The motion 
of this new plane is due to the disturbances of the planets 
and it is defined by two angles Hi and tti which depend 
on time (cf Fig 1^1). In subsection 14.31 we will determine 
accuratly the numerical expressions of Hi and tti . 

plane normal to angular momentum vector 




Fig. 2. Motion of the orbit of Venus at the date t with 
respect to the same orbit at J2000.0. 



4.1. Canonical transformations 

We denote (G' , H' , L' , g' , h' ,1') as the new set of 
Andoyer canonical variables, wich play the same role as 
(G, H, L, g,h,l) but with respect to the new moving ref- 
erence plane instead of the inertial one. We denote /' as 
the angle between the plane normal to angular momentum 
and the orbit of Venus at <. To prove that the transforma- 
tion is canonical, we have to determine a complementary 
function E depending on the new variables : 

Gdg + Hdh - Fodt = G'dg' + H'dh' - {Fq + E)dt. (8) 
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Because the tranforniation does not depend on L, I these 
two variables do not appear in ([8|) . In a given spherical tri- 
angle (A,B,C) one can show that the following differential 
relation stands (see appendix) : 

da = cos c • db + cos B • dc + sin b sin C ■ dA. (9) 

Applied to the triangle (M,N,0) (see Figl^l) this yields: 



plane normal to angular 
momentum vector 




Y oV 




plane normal to angular momentum vector 



orbit of Venus at t 



orbit of Venus at 
J2000.0 



Y oVt 



Fig. 4. Motion of the orbit of Venus at the date t with 
regard to the orbit at J2000.0. 



orbit of Venus at t 




orbit of Venus at 
N J2000.0 



Fig. 3. Details of the spherical triangle M,N,0 of FigH 
which represents the motion of the orbit of Venus at the 
date t with respect to the same orbit at J2000.0. 



dig' - g) = cos/' • d{h' - Hi) - cos J • d{h - Hi) 
-I- s\\\{h' — Hi ) sin /' • dni . 



(10) 



Using the classical equation on the same triangle : 

cos / = cos /' cosTTi — sin/' siuTTi cos(/i' — Hi). (11) 

Then, multiplying both sides by G and adding - /"o dt, we 
get : 

G-dg + H -dh-Fo-dt = G' -dg +H' -dh' -K -dt (12) 
where G = G' , H' = G' cos I' et K ^ Fq + E with 
dUi 



E = //'(I -cosTTi) 



dt 



-hG'sin/' • -sinTTi cos(/i'-ni) 



diTi 
'~dt 



sin(/i' - Hi)]. 



(13) 



(14) 



We see that the time does not appear explicitly in Fo but 
appears in the new Hamiltonian K through the variables 
TTi and Hi expressing the moving orbital plane. Thanks 
to (|12p . we have shown that the transformation from 
{G,H,L,g,h,l) to (G',H\L',g\h',l') is canonicaLThis 
canonical transformation is given in detail in Kinoshita 
(1977). To calculate the new Hamiltonian, we have to de- 
termine the values of h' — Hi, where h' is an angle con- 
structed on two different planes, the orbit of Venus at 
J2000.0 and the orbit of Venus at the date t. It is difficult 
to calculate this angle, so we can do a new transformation, 
changing h' by where hd is the sum of two angles on the 



orbit of Venus at the date t. This small transformation is 
particular to the present study and was not done in other 
works (Kinoshita 1977, Souchay et al. 1999) 

We have hd = Tovt M+MO where Tovt is the so-called 
"departure point" for Venus, which constitutes our new 
reference point. The choice of this point is justified by the 
condition of non-rotation wich characterizes it, described 
in detail by Guinot (1979) and Capitaine et al (1986). 
ToVt is the natural point to measure any motion along 
the moving plane, hd is given by the equation : 



with 



(1 — cos 7ri)c?ni 



(15) 



(16) 



To prove that this transformation is canonical, we can 
determine a function E' depending on the new variables 
such that : 



H'dh' - Kdt = H'dhd -{K + E')dt. 



(17) 



Because the tranformation does not depend on G and g, 
these two variables do not appear in ([T7]). Using equation 
([T5|) wc have: 

H'ds = -~E'dt. (18) 



Thus : 



//'^ = /?'. 

dt 



(19) 



Notice that s is a very small quantity, so that to a first 
approximation E' is negligible. Finally, we have a canon- 
ical transformation with K' as a new Hamiltonian of the 
system : 

K' ^ Fo+E + E' 

— Fq + G'sin/' • • siuTTi cos{hd — Hi) — 



dni 
'dt 



sm{hd - Hi)]. 



(20) 



4.2. Ecliptic coordinates of the orbital pole and 
reference point TOV 

To calculate the angles that characterize the moving plane 
of the orbit of Venus at the date t, a reference point on 
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the orbit of Venus at J2000.0 is needed. For the Earth the 
reference point is the vernal equinox Tq which coincides 
with the ascending node of the ecliptic of J2000.0 with 
the celestial equator of the same date. In a similar way we 
determine the coordinates of the ascending node Toy of 
the orbit of Venus at J2000.0 with respect to the equator 
of Venus at this sap e date. To calculate the coordinates 
of Toy we need Pov , the unit vector directed towards the 
orbital pole at t = defined by : 

sin io sin J7o 
Pov = I — siniocosilo 
cos io 



orbital pole of Venus 



(21) 



where io is the angle between the orbit of Venus and the 
ecliptic at J2000.0. fto is the longitude of the ascending 
node at the same date. Simon et al.(1994) give : 



io = 3.39466189 17o = 76.67992019. 



(22) 



The relative angular distance between the three poles (the 
poles of angular momentum, of figure and of rotation) is 
supposed to be very small, as is the case for the Earth. 
Throughout this study we will consider that the three 
poles coincide. So we denote P^ as the unit vector directed 
towards the pole of Venus a.t t — defined by : 



pole of Venus 




equator of Venus 



orbit of Venus 



Fig. 5. Relation between the venusian orbital pole and the 
venusian pole of angular momentum at the date t = 0. 

4.3. Determination of Hi and vri 

Using the previous section we can now determine the an- 
gles Hi and tti, which are respectively the longitude of 
the ascending node and the angle between the two orbital 
planes. The angle Hi (see Fig|4l) is equal to the sum of 
the angle v and the angle vi (see FiglHl) : 



(28) 



P, 



/cos Ip COS 6p\ 



sin Ip COS bp 



V 



( cos (5p cos 



sin Q!p cos (5p 



V 



sin(5p 



(23) 



where Zp and 6° are the ecliptic coordinates, respectively 
the longitude and the latitude of the pole of Venus. P is 
the matrix that converts equatorial coordinates to eclip- 
tic coordinates. 0:^,^^ are the equatorial coordinates of 
the Venusian north pole with respect to the equator at 
J2000.0. We take (0;°, 5^) according to the lAU report on 
cartographic coordinates(1991)(Habibullin, 1995). 




orbit of Venus 



orbit of Venus at epoch 
J2000.0 



ecliptic of J2000.0 



Fig. 6. Motion of the orbit of Venus at the date t with 
respect to the same orbit at J2000.0. 



a° = 272.76 (5° = 67.16. 



(24) 



To determine Toy we use the following equations (see 
FigEl) : 



cos /o = Po 



sin/o 



0T\ 
lOTv 



_0T\ 
|0T, 



Py A Pp. 

sin Jo 



(25) 



(26) 



where Py and Pov are the unit vectors described previ- 
ously, /q is the obliquity (angle between the axis of angu- 
lar momentum and the OZ inertial axis) and is the unit 
> 

vector along OToy. 
We finally obtain : 



Jo = 2.634, hi 



-0.103, ll = 57.75. (27) 



The relationship between tti, w, io, it, f^o and J7t are 
derived from the spherical triangle (P,L,M) : 

cos TTi — cos ig cos it -f sin ig sin it cos(r2t ^ ^^o) (29) 
coswsinTTi — — sinio cos it+cosio sinit cos(ilt^^^o) (30) 



sinusinTTi = sin it sin(f2f — f2o) 



(31) 



where io, it are respectively the inclinations of the orbit 
of Venus at J2000.0 and at the date t with respect to 
the ecliptic at J2000.0. JIq et fit are the longitude of the 
ascending node respectively at J2000.0 and at the date t. 
Simon et al. (1994) gave : 

it = 3.39466189 - 30" .84437^ - 11" .6783<2 -I- 0" .03338^^ 
VLt = 76.67992019 - 10008" .4815* - 51"3261i2 - 0".58910i^ 
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where t is counted in thousands of Julian years start- 
ing from J2000.0. In this section we choose to carry our 
calculations to third order in time. The relationships be- 
tween vi,bp and ig are derived from the spherical triangle 
(ToyGP) (see FigEl): 



orbit of Venus ai 
epoch J2000.() 




which is caused by planetary perturbations. The expres- 
sions of Fo and E + E' have been set in the previous sec- 
tion. U is the disturbing potential due to external disturb- 
ing bodies. Here the external disturbing body is the Sun 
(The perturbation due to the planets can be neglected), 
and its disturbing potential is given by : 

U = Ui + U2 



Ui = — —[ P2{smS) + — - — P2 (sin^) cos2a] 



(38) 



bit of Venus at t 



Fig. 7. Motion of the orbit of Venus at the date t with 
regard to the same orbit at J2000.0. 



n—3 



sinwi sinzo = s'mbp. (32) 
Using the equations (|^^ . (|30p . (PT|) we obtain : 
vi = 18.92 

Hi (t) = ui+v = 381348".- 901." 346 i - 0.0910822 

-0" .000759571 4- O(t^) (33) 

7ri(t) = 59".2626 i-h0.0285354 t2+o.00028575 t^+Oif^) 

(34) 

where t is counted in thousands of Julian years starting 
from J2000.0. For the Earth we have : 

ni£;(t) = 629543." - 867.927 i + 0".1534 (35) 
+0" .0000053 t^ + 0{t^) 



— ^ P™ (sin S) ■ (Cnm COS ma + Snm sin ma)] (39) 

m— 1 

where G is the gravitationnal constant, M' is the mass of 
the Sun, r is the distance between its barycenter and the 
barycenter of Venus, a and S are respectively the plan- 
etocentric longitude and latitude of the Sun (not to be 
confused with the usual equatorial coordinates), with re- 
spect to the mean equator of Venus and with respect to a 
meridian origin. The P™ are the classical Legendre func- 
tions given by : 

_ (-!)'"(! 

" ^ ' ~ 2»n! d»+™a; " ^ ' 

The Hamiltonian equations are : 



d{l',9',h') 
dK'' 



(41) 



7riE(i) =46". 9973 i-0". 0335053 
-0".00012374 t^ + 0{t^). 



(36) 



The relative motion of the orbital plane at t with respect 
to the orbital plane at t=Q is significantly greater for Venus 
than for the Earth. We used here the mean orbital ele- 
ments given by Simon et al.(1994). These elements are 
vahd over 3000 years with a relative accuracy of 10~^ and 
thus so too are our results 

5. The Hamiltonian of the system and the 
equation of the motion of rotation 

5.1. The equation of the motion of rotation 

As explained in section 21 the Hamiltonian related to the 
rotational motion of Venus is : 



In the following the prime notations used above are omit- 
ted for the sake of clarity. Using the following equations 
seen in section [3] : 

H L 
cos/=— and cosJ=— . (43) 

Gr Ct 

The canonical equations (|4T]) and (|42|) become (Kinoshita, 

(44) 
(45) 
(46) 



K" ^ Fo + E + E' + U 



(37) 



where Fo is the Hamiltonian for the free motion, E + E' is 
a component related to the motion of the orbit of Venus, 



1977) 






dl 

di ^ 


1 . ,dK 




dg_ _ 
dt 


1 , dK 


cot/—] 


dh 


1 dK 




dt ^ 


Gsin/ dl 




and 






dl 
'dt 


1 , dG 

-G^'^^'lE- 


1 dH 
sin/ dt 




I I dK 
" G^sin/ dh 


dK 

- cot I— 
dg 



(47) 
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Replacing / by J in the precedent equation yields the 
variation of J. / is the obliquity and J is the small an- 
gle between the angular momentum axis and the figure 
axis. To solve the six differential equations (pij) and 
Kinoshita (1977) used the Hori's method (1966). Here we 
solve a simplified system, i.e the potential U2 does not 
appear in our study. To know the motion of precession 
and nutation of Venus, we have to solve only the following 
equations : 



this transformation we obtain: 



dl_ 

dt 



dh 
It 



1 dH^ 



Ircot/— 
G'^'^ dt sin/ dt 



1 I dK 

G^shi7"aft" 

1 dK 
'GsmI~dT 



cot / 



dK 



(48) 
(49) 



where / characterizes the obliquity and h the motion of 
precession-nutation in longitude. These equations can also 
be written as : 



dJ 

dt 



I, dG 1 dH, 

= — cot/ 

G^ dt sin/ dt ^ 



1,1 dUi 



G' 



l/ dh 



cot / 



(50) 



dh _ 1 dUi 
H " ^Gsml~dr 



(51) 



Integrating these equations, we obtain the variations of 
the angles h and / : 



■ T»u I Uidt-cotl^ I Uidt] 
G sm I ah J og 



(52) 



Ah 



1 



d 



G sin I dl 



Uidt. 



(53) 



5.2. Development of the disturbing function 

The disturbing potential ((38)) is a functions of the modified 
Legendre polynomials P2{sm5) and P|(sini5) and of the 
planetocentric longitude a and latitude 5 of the Sun, with 
respect to the mean equator of Venus and with respect to 
a meridian origin. From the ephemerides we only know the 
longitude A and the latitude /3 of the Sun with respect to 
the orbit of Venus. So we use the transformation described 
by Kinoshita (1977) and based on the Jacobi polynomials, 
which expresses a and 5 as a function of A and fi. Applying 



P2(sin(5) =-(3cos2j-l) 



(Scos^/- l)P2(sin/3) 



--sin2/P2^(sin;3) cos2(A - h) 



- sin 2 J 



■ sin 2/P2 (sin /?) cos g 



- ^ i(l + ecos/)(-l + 2ecos/) 

e=±l 

P2 (sin/?) sin(A - h - eg) 

— ^ iesin/(l ecos/) 



e=±l 
32 



P2^(sin P) cos(2A ~2h- eg) 



+ sin^ J 



sin^ /P2 (sin /?) cos 2g + — e sin / 



f=±i 



(1 + ecos/)P2^(sin/3)sin(A -h~ 2eg) 



16 



^ (1 + e cos /)2p|(sin P) cos 2(A - /i - eg) 



e=±l 



(54) 



and 



-(3cos2/- l)P2(sin/ 



P| (sin 5) cos 2a = 3 sin^ J 

cos 2/ + - ^ sin2/P2 (sin/3) sin(A ~h~ 2el) 



e=±l 

+ \ sin^ /P|(sin P) cos 2{\ - h - el) 
8 

-|- ^ psin J(l + pcos J) 
p=±i 

3 

- - sin2/P2(sin/3) cos{2 pi + g) 

- ^ i(l + ecos/)(-l + 2ecos/) 
£=±1 ^ 

P2 (sin P) sin(A -h- 2pel - eg) 

— ^ iesin/(l + ecos/) 

e=±l 

P2^(sin P) cos(2A - 2h - 2pel - eg) 
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J2 \{l + pcosJ)' 



p=±l 



3 sin^ /P2(sin jS) cos(2/ + 2pg) 



esin/(l + ecos/) 



e=±l 



P2 (sin 13) sin(A -h- 2pel - 2eg) 
+ E J(l + ecos/)2 



e=±l 
32 



P2 (sin cos 2(A - /z, - pe^ - eg) 



Using (|^ we obtain 



P2(sin/3) = -(3sin2/3-l) 



(55) 



6. Determination of the precession and the 
nutation of Venus 

Replacing Ui in ((52|) and ((53|) we obtain : 
A/i / / r 



Gsinldl J \ 2 
1 3 

-(Scos^/- 1) - -sin2/cos2(A - /i) 
4 4 . 



A-B[3 



|^-sin^Zcos(2Z + 25) 



^ -(1+6C0S/)^ 



c=±l 



cos 2{X — h — el ~ eg) 



\dt 



(59) 



Pj(sin/3) = -3sin/3(l -sin2/3)^ 
P|(sin/3) = 3(1 -sin^/?) 

where /3 is the latitude of the disturbing body with respect 
to the orbit of Venus at the date t. In our problem the 
disturbing body is the Sun, so that (3 — 0. This yields : 



P2(sin/3) 



P2i(sin/3) = 
P|(sin/3) = 3. 

Assuming that the angle J is small (J w 0),as it the 
case for the Earth, ([M)) and (|55)) can be written : 

1 3 
P2(sin(5) = [--(3cos2/- 1) - - sin^ /cos 2(A - /i)] (56) 

and 



P2^ (sin S) cos 2a = [- sin^ / cos(2/ + 2g) 
+ ^(1 + ecos/)2 • cos2(A - h - el - eg)]. (57) 

e=±l 

Thus the disturbing potential is simplified in the following 
form : 



Ui = 



GM' 
3 



- -(3cos^/ - 1) 

4^ ' 



^ sin^ /cos2(A - h) 
A-B 



3 

- sin^ / cos(2? + 2g) 



- (1 + e cos J)^ • cos 2( A - /i - - e.g) 



e=±l 



(58) 



A/ 



1 
G 



1 d 
sin / 



GM' 



2C - ^ - P 



i(3cos2/- 1) - ^sin2/cos2(A-/i) 



yl-Pr3 



|^-sin^/cos(2/ + 25) 



+ ^ J(l+6COS/)2 
e=±l 

COS 2(A — h — el — eg) 



d f ( GM' 
cot/— ' ' 



\dt 



2C -A-B 



1 3 
- -(3cos2/- 1) - -sin2/cos2(A-/i) 
4 4 . 

A-Br3 



[-sin^/cos(2/ + 25) 



^ (1 + ecos/)^ 



e=±l 



cos 2(A — h — el — eg) 



\dt 



(60) 



To simplify the calculations, we study separately the com- 
ponents depending on the dynamical flattening and those 
depending on the triaxiality of Venus. As the dynamical 
flattening corresponds to the symmetric part at the right 
hand side of (|59l) and ([50)1 although the triaxiality cor- 
responds to the antisymmetric one, in the following the 
coefficients depending on dynamical flattening will be de- 
noted with an index "s" and the coefficients depending on 
the triaxility with an index " a" . Moreover the expression 
in dSH) and can be replaced by • where 
a is the semi-major axis of Venus given by : 



n^a^ = GM'. 



(61) 
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Finally we have the following equations 
1 d 



Ahs = 

AIS : 

and : 

Aha 



Gsinldl 

1 1 d 

G sin / dh 



Widt 



Widt 



Gsin 



1 d_ r 

J 



Wodt 



Ala^- 



where 



1 d 
sin / dh 



W2dt 



r 9 
cot/— 

dg 



W2dt 



a-^ r \ 2 

3 

— - sin^ / cos 2(A — 



-[-i(3cos^/-l) 



iX-h)]) 



(62) 
(63) 

(64) 
(65) 



(66) 



where u is the eccentric anomaly, e the eccentricity, n the 
mean motion and M the mean anomaly of Venus, we have 
the following classical development: 

lijf ^i(l + |.') + i(3e + fe>„.M 



9 



53 

-e^ cos 2M + — cos 3M. 



(72) 



For comparison we also give the corresponding coefficients 
of the development of the Earth (see Table [5]) . We have 



M 


Ls 


period 
days 


cos X 10 ^ 


f cos XlO"' 










i{l + |e^) = 5000344 [5002093] 


-96 [-105] 


1 





224,70 [365,26] 


i(3e + fE') = 101584 [250708] 


-7164 1-6305] 


2 





112,350 |182,63| 


|e= = 1032 [6282] 


-146 1-316] 


3 





74.900 [121.753] 


f e-" = 200 [308] 


-40 [-23] 



Table 2. Development of 5 (7) of Venus (the corre- 
sponding values for the Earth are in square brackets), t 
is counted in Julian centuries. 



W2 



sin^ /cos(2/ + 2g) 



GM' a^/A- B 
a3 V 4 L2 

£=±1 

cos2(A — h ~ el — eg) 



We adopt the following notations : 

Ai/; = A/i = {Ahs + Aha) 
Ae = A/= (A/s + A/a). 



(67) 



(69) 



As the rotation of Venus is retrograd the convention used 
is not the same as lAU to the Earth. 



7. Development of the disturbing function 

Considering the level of accuracy and the time interval 
involved (3000 y), truncating the function e to first order 
in time is a sufficient approximation. Therefore we take 
the eccentricity of Venus from Simon et al.(1994) as : 



e = 0.0067719164- 0.0004776521t. 



(70) 



7.1. Terms depending on the dynamical flattening 



To solve the equations (15^ and which give the nuta- 
tion coming from the dynamical flattening, it is necessary 
to develop 5(7)^ and (7)^cos2(A - h) with respect to 
time (i.e with respect to the mean anomaly and the mean 
longitude of Sun). Using Kepler's equation : 



the following definition : 

2A = 2(07+1;) = 2{lJ+M+{v-M)) = 2Ls+2{v-M) (73) 

where A is the true longitude of the Sun, ZJ the longitude 
of the periapse of Venus, M the mean anomaly and L5 
the mean longitude of the Sun. From Kepler's equation 
we have the following classical development : 

v-M = (2e- — )sinA/-f(-e^ e'^)sin2AfH e^sin3Af. 

^4^ ^4 24 ^ 12 

(74) 

From the trigonometric equation : 

cos2A = cos2i5Cos2(w-Af)-sin2Lssin2(t;-M) (75) 
we obtain : 

(-)3cos2(A-/i) = (-)3(cos2i:scos2(t>-M) 
r r 

-sm2Lssm2{v ~ M)). (76) 

Using the developments ([7^ and (1751) and the classical 
trigonometric relationships we obtain a development of 
(f)3cos2(A-/i) : 

(-)3cos2(A-/i) (1 - -e2)cos2Ls 
r 2 

7 123 

+ (oe-^e3)cos(2L5 + M) 
z Id 

+ (-^e+^e3)cos(2Ls-M) 
+ {^e^)cos{2Ls + 2M) 



845 
'1^ 



cos(2Ls + 3M) 



-^e^ cos(2Ls - 3M). 



u — e sin u = n{t — to) ~ M 



(71) 



(77) 
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The coefficients of tliis development are given in Table 
[3l Tfie periods in Tables [2] and [3] have been computed 



{-Y cos 2(A - /i + $) (1 - -e^) cos(2is + 2$) 



M 


is 


period 
days 


cos X 10 ' 


tcos xlO"' 





2 


112,35 [182,62] 


(1 - 1 e^) = 9998853 [9993025] 


161 [351] 


-1 


2 


224,70 [365,22] 


(-ie + ie^) = -33859 [-83540] 


2388 [2101] 


1 


2 


74,90 [121,75] 


(|e - J^e-') = 236993 [584444] 


-16718 [- 14713] 


2 


2 


56,17 [91,31] 


(^e^) = 3898 [23730] 


-550 [-1193] 


3 


2 


44.94 [73,05] 


Sfe3=54 [821] 


-4 [-20] 


-3 


2 


-224.70 [-365.25] 


= 0[0] 


0[0] 



Table 3. Development of (^) cos(2(A - h) of Venus (the 
corresponding vahies for the Earth are in square brackets) . 
t is counted in Julian centuries. 



through the mean elements of the planets given by the 
ephemerides VSOP87 of Simon et al.(1994). Because the 
eccentricity of Venus is smaller than the eccentricity of 
the Earth, the coefficients of the development for Venus 
are smaller than the corresponding ones for the Earth, 
except for the leading one in Table [H 



7.2. Terms depending on the triaxiality 

To solve the equations ([M]) and ([55]) it is necessary 
to develop i(f)3cos2$, (f)3cos(2(A - h) - 2$) and 
(-)^ cos(2(A — h) + 2$) with respect to the mean anomaly 
M, the mean longitude Ls of Sun and to the angle 
I + g ^ ^ according to ([6|). Using classical trigonomet- 
ric equations and developments similar to those used in 
section FTTl we obtain the following : 



4 



123 



16 

1 1 



e^) cos(2Ls + M + 24>) 



H'^e + —e') cos(2Ls - M + 2$) 

+ (ye2) cos(2is + 2M + 2$) 

+ ^e^ cos{2Ls + 3M + 2$) 
48 

+ ^e^ cos(2Ls - 3M + 2$) 



1(^)3 cos(2$) 



(79) 



= (i + ^e2)cos(2<i>) 

Q 27 
+ (4e-)--e3)cos(Af-2$) 

3 27 

+ ^4''"^ 32^'^^°'^^'^ + ^'^^ 

+ le^ cos(2M - 2$) 
8 

+ de^)cos{2M + 2<^>) 
8 

53 

+—e^ cos(3M-t-2$) 
16 

53 

+— cos(3M-2$). 
16 



These are given in Tables I4I5I and [S] 



(80) 



M 


Ls 




period 
days 


cos X 10 ^ 


t COS X 10 ^ 








2 


-121.51 


(i + |e^) = 5000344 


-48 


1 





2 


-264.6 


(|e + fe'') = 50792 


-3582 


1 





-2 


78.86 


+ if/') = 50792 


-3582 


2 





2 


1490.35 


ih') = 516 


-72 


2 





- 2 


58.37 


(fe^) = 516 


-72 


3 





2 


195.26 


ffe^ = 10 





3 





-2 


46.34 


^e'* = 10 

16 






'cos2(A-/i-$) (1 - -e^)cos(2Ls - 2$) 



123 
16 



e^) cos(2L5 + M - 2$) 



Table 4. Development of 
counted in Julian centuries. 



+ (-^^^ + 4^') cos(2L5 -M- 2$) 
2 Id 

-^(ye^) cos(2L5 + 2M - 2$) 

+— cos(2i5 + 3Af - 2$) 
48 

-F-J-e^ cos(2is - 3Af - 2$) 
48 

(78) 



' cos(2<i>) of Venus, t is 
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M 


Ls 


$ 


psriod 
days 


cos X 10 " 


t cos X 10 " 





2 


- 2 


58.37 


(1 - |e^) = 9998853 


161 


-1 


2 


- 2 


78.87 


(-ie+ ie=) = -33859 


2388 


1 


2 


- 2 


46.34 


{|e - ^e-'') = 236993 


-16718 


2 


2 


- 2 


38.41 




-550 


-3 


2 


-2 


264.59 


= 





3 


2 


-2 


32.80 


l^e^ = 54 


-4 



Table 5. Development of (f )^ cos(2(A-ft.)-2$) of Venus. 
t is counted in Julian centuries. 
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M 


Ls 




days 


cos X 10 ^ 


tcos xlO"^ 





2 


2 


1490.35 


(1 - |e^) = 9998853 


161 


-1 


2 


2 


-2C4.6 


{-ie-F -Le') = -33859 


2388 


1 


2 


2 


195.26 


( Je - If = 236993 


-16718 


2 


2 


2 


104.47 


(^e^) = 3898 


550 


-3 


2 


2 


224.70 







3 


2 


2 


71.31 




-4 



Table 6. Development of (f )^ cos(2(A-/i)+2$) of Venus. 
t is counted in Julian centuries. 



We note that the coefficients in Tables [S] and [S] are the 
same as the coefficients of Table [31 because the calcula- 
tions are similar and they do not depend on the frequency 
associated with each argument, but only on the eccentric- 
ity of Venus. Nevertheless the corresponding periods are 
different because their calculation includes the argument 
<i>, the sideral rotation of Venus given by : 

2TT 

'^^l + g^-ijr) t (81) 

where Ty is the sidereal period of the rotation of Venus. 
The numerical value for Ty is Ty = 243.02 days. Notice 
also that because the rotation is retrograde, $ has a nega- 
tive sign, and that the period of Ls (225 d) is close to 
the period of $ (243 d). Therefore the period for the 
leading componant 2Ls + 2$ in Table M (1490.35 d) is 
much longer than the period of the leading component 
2Ls - 2$ in Table [5] (58.37 d). We have Ls = ujo t, 
M = Uouj t and $ = t where Wo = (2n/224.70080)rd/d, 
ujou = (2n/224.70082)rd/d and Ur = -(2n/243.02)rd/d. 
The very small difference between and ljJou is due to 
the longitude of periapse Zo which corresponds to a very 
low frequency. As a consequence we do not change our 
calculations significantly by considering that Ls — M 
{Ls = M + uj). 

8. Precession and nutation which depend on the 
longitude and obliquity 

In the following, we calculate the precession and the nuta- 
tion of Venus depending on both the dynamical flattening 
and the dynamical triaxiality. 



8.1. Precession and the nutation of Venus depending 
on the dynamical flattening 

According to and we compute the precession- 
nutation in longitude (Ahg) and obliquity (A/^) due to 
the dynamical flattening : 

1 d f GM' \ f2C - A- B\— ^ 



sin / dl 
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AIs = 



1 d f GM' 1 f2C-A-B 



sin / dh 
d 

sin / dh 



Widt (83) 



Widt 



where 



Wi ^ - —{3 cos^ / - 1) - - sin^ / cos 2(A - h)]{8A) 



12 

Ahs and A/^ in 
scaling factor Kg : 



and 



depend directly on the 



GM' 



Ga3 



2C -A-B 



Using the third Kepler 's law we obtain : 
,2C - A - B, ir? 



2C 



(85) 



(86) 



where n is the mean motion of Venus. The dynamical flat- 
tening, defined as By — , is a dimensionless pa- 
rameter. We supposed here to the first order of approx- 
imation that the components wi and W2 of the rotation 
of Venus are negligible with respect to the component W3 
along the figure axis (0,z), as is the case for the Earth 
(Kinoshita, 1977). So we have 



G = Auji + Buj2 + Ccjs w Guj 



(87) 



with the values of A,B,G adopted in section [TU] and that of 
the dynamical flattening given which is deduced directly 
from them (given in Table [1} , we find : 



Ks = -8957".55± 133. 29/ Ju/mn. cy. 



In the case of the Earth, the scaling factor due to the 
gravitational influence of the Sun is (Souchay et al.l999) 



= 3475". 36/ Julian, cy. 
Then we have the following result : 



Ks _ ,Tf_. ,i£E_. f Hy . 



-2.577. 



(89) 



(90) 



Thus for Venus, the value of the constant of the precession- 
nutation due to the Sun is roughly 2.6 times bigger that 
the corresponding one for the Earth. Through the devel- 
opment done previously (cf section[7]) and according to the 
conventional notations and we get the preces- 

sion and the nutation of Venus depending on its dynami- 
cal flattening. We will see in the following sub-section that 
the triaxiality does not contribute to the precession. So we 
obtain the following result for the precession : 



iP^Ks cos I j {l + -e^)dt = 4474" .35t - 0.021t^ (91) 

where t is counted in Julian centuries . The nutation coef- 
ficients in longitude and in obliquity are given respectively 
in Tables [7] and [S] All our results are given with respect 
to the moving orbit at the date t. 



Table 7. A^Ps = A/i^: nutation coefficients in longitude 
of Venus depending on its dynamical flattening. 



Argument 


Period 
days 


sin 
arc second 
(10-7) 


t sin 

arc second/ Julian cy 
(10-7) 


cos 
arc second 
(10-7) 


2L, 


112.35 


21900468 


352 





M 


224.70 


-889997 


62765 


61 


2Ls + M 


74.90 


346057 


- 2 4412 


-7 


2L, - M 


224.70 


-148323 


10461 


10 


2L, + M 


56.17 


4269 


- 602 





2M 


112.35 


-4521 


640 






Table 8. Ae<; = A/s : nutation coefficients in obliquity of 
Venus depending on its dynamical flattening. 



Argument 


Period 
days 


cos 
arc second 
(10-7) 


t cos 

arc second/.! ulian century 
(10-7) 


sin 
arc second 
(10-7) 




112.35 


-100741 


-16 





2L, + M 


74.90 


-15919 


1123 





•iLs - M 


224.70 


6822 


-481 





2L, + 2M 


.56.17 


-196 


27 






These coefficients are of the same order of magnitude 
as those for the Earth due to the action of the Sun found 
by Kinoshita (1977) and Souchay et al. (1999). The largest 
nutation coefficient has roughly a 2" amplitude whereas 
for the Earth it has a 1" amplitude. As in the case of the 
precession, this difference comes from the larger value of 
the constant Kg and is slightly compensated for after inte- 
gration by the fact that the period of revolution of Venus 
around the Sun (225 d) is shorter than that of the Earth. 
The largest nutation coefficient in obliquity is compara- 
tively very small with a 0.1" amplitude whereas for the 
Earth it has a 0.5" amplitude. The difference is due to the 
small obliquity of Venus (/ w 3), our results depending on 
sin/. 

8.2. Precession and nutation of Venus depending on 
dynamical triaxiality 



Starting from (|64| and (j65l) , we can compute the nutations 
Afia and Ala respectively in longitude and in obliquity, 
due to the triaxiality of Venus : 



Aha = - 



1 d 
sin / dl 

Kg d 

sin / dl 



f GM' 1 fA-B\—^ 



(92) 



W^dt 
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Ala 



1 d 
sin / dh 



Kg d 

sin / dh 



GM' 1 / A- B 



Wodt 



GM' 1 A-B—^ 
3 — — — : — W2dt 



G 



W2dt - K a cot I 



dg 



Widt 



(93) 
(94) 



where 

^ (-)3 [1 sin^ / . cos(2$) 
r '-2 

+ ^ -{l + ecosif 

£=±1 

•cos2(A -h- e$)]. 



(95) 



A/ia and A/a in and ((M)) depend directly on the 
scaling factor Ka : 



ifa = 3- 



.GM' 



(96) 



where M', G, G and a have been defined previously. Still 
using Kepler's third law we get : 



A-B, 



(97) 



The triaxiality, defined as Ty — , is also, as Hv^ a 
dimensionless parameter, here we assume again that : 



G« Cw. 



(98) 



Then, with the value of the triaxiality (given in Table [T]), 
we get : 



Ka = 1131".23± 17.09 Julian/ cy. 



(99) 



In the case of the Earth, the scaling factor due to the 
gravitational infiuence of the Sun is (Souchay et al.l999) 



iff = -5" 68 Julian/ cy. 
Thus we have the following ratios : 



|i-(4)-(^)-(^) 

nru LU 1e 



-199.16. 



(100) 



(101) 



This result shows that the coefficients of nutation due to 
the triaxility are considerably larger for Venus than for 
the Earth. Thanks to the development done previously (cf 
section HTTj) and according to the equations (|^ and 
we get the nutation of Venus depending on its triaxiality. 
The nutation coefficients in longitude and in obliquity are 
given respectively in Tables IH] and [TUl All our results are 
given with respect to the moving orbit at the date t. 

For Venus, the largest nutation coefficient in longitude 
for the terms depending on its triaxiality has a 0" .6 ampli- 
tude and the largest nutation coefficient in obliquity has 
a 0" .03 amplitude. Their period corresponds to half the 
period of rotation of the planets, 121.5 d. 



Table 9. A^Pq = Aha- nutation coefficient in longitude of 
Venus depending on its triaxiality. 



Period sin t sin cos 

days arc second arc second/Julian cy arc second 







10^7 


10^7 


10^7 


2* 


-121.51 


-5994459 


58 





2Ls - 2$ 


58.37 


-2880826 


-46 





M + 24' 


-264.6 


-132590 


9351 


- 11 


A/ - 24> + 2Ls 


46.34 


-54201 


3823 





M -2i 


78.86 


39519 


- 2787 





2Ls + 2* 


1490.35 


38866 








-M - 2* + 2Ls 


78.86 


13179 


-929 





2M + 24> 


1490.35 


7587 


- 1058 


-7 


2M - 2* + 2Ls 


38.41 


-739 


104 





2M - 24> 


-58.37 


297 


-41 





A/ + 24> + 2Ls 


195.26 


121 


- 8 





-Al + 24 + 2Ls 


-264.66 


23 


-2 





2M + 2* + 2Ls 


104.47 


1 









Table 10. Aea = A/q mutation coefficients in obliquity 
of Venus depending on its triaxiality. 



Argument 


Period 
days 


cos 
arc second 


t cos 

ire second / Julian cy 


sin 
ire second 


24> 


-121.51 


275453 


-3 





2ls - 24 


-58.37 


-132365 


-2 





A/ + 24 


-264.6 


6093 


-430 





A/ - 24 4- 2Ls 


46.34 


-2491 


176 





21s + 24 


1490.35 


-1786 








A/ - 24 


78.86 


-1816 


128 





-a; - 24 + 2Ls 


78.86 


606 


-43 





2M + 24 


1490.35 


-348 


47 





2A/ - 24 + 2Ls 


38.41 


-35 


5 





2M - 24 


-58.37 


-14 


2 





A/ + 24 4- 2Ls 


195.26 


-6 








-M + 24 + 2Ls 


-264.6 


-1 








2M + 24 + 2Ls 


104.47 












From the Tables [7] to [TUl we choose to calculate nu- 
merically from their analytical expression the precession- 
nutation with respect to i, over a period of 4000 d. Fig|Sl 
represents the combined precession and nutation of Venus 
in longitude. The periodic part stands for the nutation and 
the linear part for the precession (see ([91])). In Fig[9| and 
Fig llOI . we represent the nutation of Venus respectively 
in longitude and in obliquity depending on its dynami- 
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Fig. 8. The precession and the nutation of Venus in lon- 
gitude for a 4000 day time span. 
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Fig. 9. The nutation in longitude of Venus depending 
on its dynamical flattening (cyan) and depending on its 
triaxiality (blue) for a 4000 day time span, from J2000.0. 
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Fig. 10. The nutation in obliquity of Venus depending 
on its dynamical flattening (cyan) and depending on its 
triaxiality (blue) for a 4000 day time span, from J2000.0. 

cal flattening (black curve) and its triaxiality (red curve) . 
In FigfTTl and FigH^l, we represent the nutation respec- 
tively in longitude and in obliquity when combining the 
two components above. In Figsl8]to[T2l we choose the 4000 
d time span to clearly see the leading oscillations. Recall 
that our study is valid over 3000 years. The nutation of 
Venus depending on its triaxiality is dominated by a si- 
nusoid with a period of 121.51 d and another one with 
a period of 58.37 d. On the other hand the nutation de- 
pending on its dynamical flattening is largely dominated 
by a single sinusoid with a period of 112.35 d. We can also 
remark (FiglHl and FigfTUl) that the nutation coefficients 
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Fig. 11. The nutation of Venus in longitude for a 4000 
day time span, from J2000.0 
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Fig. 12. The nutation of Venus in obliquity for a 4000 
day time span, from J2000.0. 

of Venus depending on its triaxiality are of the same or- 
der of magnitude as the nutation coefficients depending on 
its dynamical flattening, whereas for the Earth they are 
small and negligible in comparison ( Kinoshita,1997 and 
Souchay et al. 1999). Two reasons explain this : 

— the ratio of scaling factors is ^ = « ^ for Venus, 
whereas for the Earth it is two orders of magnitude 
smaller : = « ^ 

— the frequency of the sidereal angle $ ss / + 5 which 
enters in the denominator during the integrations in 
([92]) and ((93|) is roughly 243 times larger for the Earth 
than for Venus 

Thus we have pointed out a significant difference between 
the nutations of the two planets. 

9. Motion of the pole of Venus in space 

From the previous section, it is possible to describe the 
motion of the pole of Venus with respect to its or- 
bital plane at the date t. The rectangular coordinates 
{Xp, Yp,Zp) of the pole are given by : 

Xp\ /- sin(/o + A/) sin(/i + A/i) \ 
Vp = - sin(/o + A/) cos(/i + Ah) (102) 
ZpJ \ cos(/o+A/) / 

The precession in longitude ^ is given by 

= *t = -/i = 4474" .35i ±66.5 {t in Julian/ cy) (103) 
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Iq is the nominal value of the obliquity 



with these values and adopting 



c 
WW 



— 0.3360, we get: 



lo = -eo = -2.63. 



(104) 



Replacing A/ and Ah by their numerical values from 
Tables [71 to [TUl we show in FigH^l the motion of the Venus 
figure axis in space over a one century time interval, for 
which the combined motion of precession and nutation 
appears clearly. 

From (|103p we directly determine the period of preces- 
sion of Venus, i.e 28965. 10 ±437 years. It is slightly longer 
than the period of precession of the Earth, i.e. 25712.4 
y (Bretagnon et al. 1997). Notice that in the case of the 
Earth, the Moon contributes roughly 2 /3 and the Sun 1 /3 
to the precession rate. In the case of Venus, only the Sun 
contributes significantly to the precession rate. However, 
as we have shown in section [51 this solar contribution is 
roughly 2.5 times greater than the corresponding one for 
the Earth. This explains why the periods of precession are 
rather equivalent. 




Fig. 13. (X,Y) motion of Venus polar axis in space for a 
one century time span 



10. Numerical values of 
comparison of results 



2C-A-B 



2C 



and 



A-B 
AC 



and 



From the various calculations done in the previous sec- 
tions, we know that the accuracy in the determination of 
the precession rate and of the coefficients of nutation 
depends directly on the quality of determination of the 
dynamical ellipticity Hy = ^'~^~2c~^ ^'^'^ triaxiality 
Tv = Therefore for this purpose it is crucial to ob- 

tain accurate values of the moments of inertia. The value 
of ratio -j^^ is inperfectly known. Kozlovskaya (1966) 
obtained a value in the [0.321;0.360] interval, Shen and 
Zhang (1988) in the [0.321;0.350] one and Yoder(1995) in 
the [0.331;0.341] one. Williams (private communication) 
gives the following values relative to the differences of the 
moments of inertia: 



C-A 

C 

B-A 



C 



= 1.643 X 10~^■ 
6.631 X 10"^. 



C-B 



C 



= 9.79 X 10-^ 



(106) 



All our calculations in this paper have been done with 
these values. When studying the rotation of Venus, 
Habbibulin (1995) has taken jj^ = 0.340 and the val- 
ues of the Stokes parameters obtained by Williams (1983) 



A ^ r> 

= 2.3324 X 10"^ — — = 1.2618 x 10"^ 



C 

B-A , 
= 1.0706 X 10"^ 



c 



(107) 



In Table [TT] we compare our results with those ob- 
tained by directly using the values of Habbibulin (1995) 
for the physical parameters of Venus, thus showing large 
differences. 

Table 11. Comparison between our results computed 
with the recent values of the moment of inertia and those 
obtained with the values of Habbibulin et al.(1995) 



Principal parameter 


Habbibulin(1995) 


This study 


Triaxiality : 


-2.67* lO""* 


-1.66* 10"'' 


Dyn. flattening : 


1.70* 10"' 


1.31 * 10"'' 


kY 


-12275". 32/cj/ 


-8957".55/ci/± 133.29 


kX 


1829".34/ci/ 


1133".28/ci/ ± 17.09 


Precession 


6133". 77i 


4474". 35t ±66.5 


Period of precession 


21156 years 


28965.10 ±436.99 years 


Largest nutation coeflacient 


3". 00 


2". 19 



Table [TTl shows that the values of the moments of iner- 
tia are crucial for our computations. The comparison with 
the final results computed by Habbibulin (1995) is not pos- 
sible because in his paper different angles, parametrisation 
and more complex developments were used. In another 
study, Zhang (1988) found a precession rate of 7828". 6/cy 
corresponding to a period of 16555 years. Altough no pre- 
cise information could be extracted from this paper, the 
difference is probably due to different values of the mo- 
ments of inertia. What we highlight here is that different 
values of moment of inertia yield fairly different results for 
the precession and the nutation as shown in Table 1111 



C-A 
B-A 



= 5.519 X 10 



= 2.228 X 10' 



-6 



C-B 
MK^ 



= 3.290 X 10" 



(105) 



11. Conclusion and Prospects 

In this paper we investigated by the rotation of Venus. 
We adopted a well suited theory of the rotation of a rigid 
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body set up by Kinoshita (1977) and developed in full de- 
tail by Souchay et al.(1999) for the Earth. In a first step 
we defined the reference frames to the various polar axes 
of the planet (axis of angular momentum, rotation axis, 
figure axis) and on its moving orbital plane. We have also 
given in full detail the parametrization of Venus rotation 
starting from the set of Andoyer canonical variables with 
respect to the moving orbital plane. In particular we have 
determined precisely the motion of the orbit of Venus at 
t with respect to the orbit at J2000.0, through the pa- 
rameters Hi and tti, for which we have given polynomial 
expressions. Moreover we have checked the value of the 
Venus obliquity (263). 

Then we calculated the disturbing function due to 
the gravitational interaction with the Sun. Applying 
Kinoshita's Hamiltonian analytical developments, we cal- 
culated the precession constant of Venus, with a precision 
and an accuracy better than in previous works (Habillulin, 
1995, Zhang, 1988). Our value for the precession in longi- 
tude is tp = 4474" .35<:/cy ± 66.5. which is more than two 
times larger than the corresponding term for the Earth 
due to the Sun {ip — 1583".99/cy) and slightly smaller 
than the combined effect of the Moon and of the Sun for 
the Earth {ip = 5000.3" /cy). We have shown that the ef- 
fect of the very small value of the dynamical ellipticity of 
Venus {Hv = 1.31 x 10~5), which should directly lower 
the amplitude of its precession and nutation, is more than 
fully compensated for by its very slow retrograde rota- 
tion. Moreover, one of the specificities of Venus is that 
its triaxiality {Ty = —1.66 x 10~6.) is of the same or- 
der as its dynamical ellipticity (see value above), unlike 
what happens for the Earth for which it is considerably 
smaller. Consequently, we have performed a full calcula- 
tion of the coefficients of nutation of Venus due to the 
Sun and presented the complete tables of nutation in lon- 
gitude (A^*) and obhquity (Ae) due both to the dynamical 
ellipticity and the triaxiality. These tables have never been 
presented in previous works. We think that this work can 
be a starting point for further studies dealing with Venus 
rotation, for which it has set up the theoretical founda- 
tion (parametrization, equations of motion etc.), such as 
a study of precession-nutation over a long time scale, the 
calculation of Oppolzer terms, the effects of the atmo- 
sphere. 



we obtain: 

sin a • da = (— sin b cos c + sin c cos b cos A) • db 
-|-(— cos b sin c + sin b cos c cos A) • dc 
-f (— sinb sin c sin A) • dA. 



(108) 



Dividing by - sin a we have the following equation: 

/ sinbcosc — sin c cos b cos A \ 
da = ■ • db 



sm a 

cos b sin c — sin b cos c cos A 



sm a 
sin b sin c sin A ^ 



dc 



sm a 



■ dA. 



(109) 



In the spherical triangle A,B,C the following equations 
holds : 



sin c sin a sin b sin c sin A 



sin C sin A 



= sin b sin C. 



Thus 



sin a cos C = cos c sin b — sinb cose cos A 

sin b cos c — sin c cos b cos A 



cos C 



sin a cos B 



sm a 



cos b sin c — sin c cos c cos A 
sin c cos b — sin b cos c cos A 



cosB 



sm a 

Replacing in (|109p . wc obtain ([9]). 
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